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ABSTRACT 

We present the first Spitzer-IRS spectral maps of the Herbig-Haro flow GGD 37 detected in lines 
of [Ne III] [O IV], [Ar III], and [Ne V]. The detection of extended [O IV] (55 eV) and some extended 
emission in [Ne V] (97 eV) indicates a shock temperature in excess of 100,000 K, in agreement with 
X-ray observations, and a shock speed in excess of 200 km s _1 . The presence of an extended pho- 
toionization or collisional ionization region indicates that GGD 37 is a highly unusual protostellar 
outflow. 

Subject headings: ISM: jets and outflows. ISM: Herbig-Haro objects. ISM: clouds 



1. INTRODUCTION 

Protostellar cores of all masses eject material in the 
form of detectable bipolar outflows when their accre- 
tion rates are sufficiently high, ie. > 10 _T M yr _1 . 
Very massive stars may trigger additional star formation 
through these powerful flows. Even the much more nu- 
merous outflows from low mass stars may impact their 
local environment with high velocity flows and stir up 
turbulence. Intermediate mass stars are an interesting 
middle scenario: are they numerous enough and powerful 
enough to trigger further star formation, or suppress it? 
One way to determine their effect on their environment 
is to analyze the physical conditions of the surrounding 
material before and after it is impacted by shocks tied to 
t he j et / outflow . 

Sho cks were categorized by their physics in IDraind 
(| 19801 ): J-type shocks in which the magnetic field from 
the shocked ions is either frozen into the shock, or is non- 
existent; J-type shocks with magnetic precursors (or ra- 
diative precursors), in which the neutral fluid undergoes 
a discontinuous "jump" in density and temperature; and 
C-type shocks in which the fluid changes continuously 
over the shock boundary. The conditions in these shocks 
vary and can be constrained through observations of fine 
structure line emission with differing appearance poten- 
tial and critical density. At 700 pc distant, the Cepheus 
A cloud complex contains a number of Herbig-Haro ob- 
jects - natural laboratories for shock models - housed 
within the Cep OB3 association. 

Cep A contains a collection of luminous protostellar 
sources; additionally these protostars are associated with 
at least two broa dly extended CO outflows, C ep A East 
and Cep A West ([Hughes fc Wouterlootll 19821): the latter 
some times referred to as GGD 37 ( Gvulbudagh ian et 
[19781 ) or HH 168 (jReipurth fc Rag a 1999) . Both flows 
have been mapped by Spi tzer-IRS ([Sonnentrucker et"aTl 
120061 : iNeufeld et aLl l2006a). The e xtinction toward GGD 
37 i s generally thought t o be small ( Fro ebrich et al.l l 2002) 
but [Wright et al. (1996) apply an Ay of 16 to their data. 
High-resolution near-IR imaging ([Raines et al. 1 120001 ) of 
the molecular emission from GGD 37 reveals a chain of 



bow shocks in GGD37. The SE-NW orientation of most 
of these shocks suggest an outflow driven from the vicin- 
ity of HW3c (a 10 M protostar in Cep A East), while 
other bow shocks suggest a flow orientation that is more 
east- west. Additionally a chain of several extremely high 
proper motion (~ 800 km s _1 ) knots of [Fe II] 1.644 /im 
emission connects the SE-NW flow but is notably ab- 
sent in the E-W flow. The [Fe II] emission peaks around 
a trio of radio sources, centered at "W2" ([Gar ay et al.l 
1996) (see Figure [1]) . 

The Cep A region has been m apped in both hard and 
soft X-rays w ith XMM-Newton (Pravdo & Tsuboi 20051 : 
iPravdo et aTll2009[ ) and Chandra (Schneid eret al.ll2009[ ). 
Both HW2 and HW3c are detected in hard X-rays, while 
GGD 37 (incl uding the W1-W2-W3 re gion) is detected 
in soft X-rays. IPravdo fc Tsuboil (j2005h note several dis- 
tinct soft X-ray sources within the GGD 37 flow, one near 
the W2 region, and one at the far end of the flow. They 
conclude for GGD 37 that L x = 3 x 10 30 erg s" 1 , and 
infer a shock velocity of 620 km s _1 , and a temperature 
of 4 x 1Q 6 K from comparisons to L1551-IRS 5 / HH 
154 (jFavata et aLll2QQ2t iBaHv et al.ll2003D . They observe 
that the flow is extreme in temperature, although not 
in luminosity. The X-ray emissions are offset from both 
the optical and radio emissions, and they conclude this is 
due to the complex morphology of the region and multi- 
ple driving sources, as well as anisotropy in the medium. 
In this paper we present Spitzer-IRS maps of highly ion- 
ized mid-IR forbidden lines in GGD 37 and show that 
the offset and extended nature of these Spitzer lines, and 
the X-ray emission, are directly attributable to the cloud 
shock. 

2. OBSERVATIONS 

2.1. Data Reduction 

We utilized the InfraRed Spectrograph (IRS) onboard 
Spitzer to observe GGD 37 as part of our guaranteed and 
general observing time (program 2, PI: J.R. Houck; pro- 
gram 30167, PI: D. A. Neufeld); for the full dataset, see 
Green et al. (2010a, in prep.). We observed this region in 
a Spitzer-IRS map grid consisting of two adjacent rectan- 
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gles, on two separate occasions; first in Short-High (SH; 
spectral resolution R ~ 600) and Long-High (LH; R ~ 
600) in 2004, and again in Short-Low (SL; R - 60-120) 
and in deeper exposures in Long-High in 2007. These 
maps cover approximately the same regions of the sky, 
although Cep A East was observed concurrently with 
the first set of observations of GGD 37. The resolution 
of individual pixels of the IRS is given in the Spitzer Ob- 
servers' Manual Q as ranging from ~ 2 - 5 arcsec squares, 
depending upon the module in question; therefore the 
maps contain differing spatial resolution, although they 
have all been regridded and oversampled for display pur- 
poses. 

For full details of the data reduction process, see Green 
et al. (2010a); here we briefly summarize the procedure. 
The basic calibrated data (BCD) were processed at the 
Spitzer Science Center (SSC) using version S12.0, S15.3, 
or SI 7. 2 of the processing; pipeline , and then reduced 
using SMART (jffigdon et all l2004f ). modified by addi- 
tional routines that we have developed to deal with map 
grids (as opposed to single observation staring mode). 
We process the data beyond the traditional SMART re- 
duction by removing bad pixels in all modules, using a 
"grand" rogue mask created from a superposition of the 
bad pixel masks of each separate observing campaign. 
Rogues that occur once at the 4<r level are considered 
to be permanently bad pixels and thus the grand rogue 
mask contains 25% more rogue pixels than a single cam- 
paign rogue mask. Bad pixels are fixed using a nearest 
neighbor routine (imclean.pro and imnan.pro). 

Next we extracted individual spectra not simply from 
each slit position, but from each individual resolution el- 
ement of each slit position. In the case of SH and LH, we 
produce five (semi-)independent spectra, spatially sepa- 
rated by the resolution of a single pixel in the array; in SL 
we use 32 independent positions to generate spectra. The 
spectra are then re-gridded onto a regular grid in right 
ascension and declination. Pixels with (nearly) identical 
positions are averaged together. We apply the slit loss 
correction function (J. D. Smith, private communication) 
to correct for flux- calibration created for point-sources, 
used by the SSC. The result is a downward correction to 
the flux to ~ 60-90% of the original value. 

2.2. Extended Detection of Lines with Appearance 
Potential > 13.6 eV 

We report the first detection of [Ar III], [Ne III], [O 
IV] or [Ne V] in a Herbig-Haro flow. The [O IV] (25.89 
jam) detection is fit separately both spectrally and spa- 
tially from the much brig hter nearby [Fe II] line (25.99 
/am). iGuiles et al.l ([2007) observed emission from ap- 
parent spectral lines several resolution elements on both 
sides of extremely bright lines in their sources, but they 
were ascribed to data artifacts due to the symmetric na- 
ture of the "lines" around the central peak. In this case 
the [Fe II] line brightness is not extreme and there is no 
corresponding bump on the long-wavelength side of the 
[Fe II] 26 /am line. The [Ar III] map is co-spatial with 
the other lines, with the exception that it appears to con- 
tribute to the westernmost clump in the map, unlike the 
other species of higher appearance potential. Addition- 
ally there is a possible detection of [Ar III] (8.99 /im) 

1 http:/ /ssc. spitzer. caltech.edu/documents/som/ 



but it is suspect due to lack of coherence in its mapped 
spatial structure (Green et al. 2010b, in prep.). 

The [Ne V] emission is confined to the area around W2 
and a point at the northwest end of the flow (Figure [2j), 
roughly coinciding with t he edge o f the lower excitation 
line emission from [Fe II]. iRainesI (j2000l . his figure 3.15) 
noted that the peak of the [Fe II] (1.644 /im) emission 
was shifted to the east by ~ 3" from the peak of the H2 
(2.122 jam) emission; we note a similar result in our maps 
(Figure [3]) - although the spatial resolution is consider- 
ably lower than in Raines' images - and we can clearly 
see that the peak emission from all of the ionized species 
is located east of the peak of the molecular emission. 
Additionally, the higher ionization species appear to be 
themselves shifted eastward by a further ~ 5" compared 
to the lower ionization species, suggesting that we are 
resolving the postshock region. 

Many other species, includ ing [Ne III, [ Si Hi, and [S I] 
were detected in GGD37 in Wrigh t et all (1996) with the 
Infrared Space Observatory (ISO); however strong lines 
in the IRS data such as [Ne III] were not noted in the 
ISO data; this may be due to the large beam size of ISO 
and the small spatial extent of the [Ne III] emission. 

3. DISCUSSION 

Originally detecte d as a radio emitter by 
iHughes fc Wouterlootl ([1982) and resolved int o three 
separa te components (Wl, W2, and W3) by iHughesI 
(jl989[ ) using VLA data, W2 is a prime candidate for 
an intermediate mass protost ellar object drivi ng an 
exceptionally powerful outflow (jGarav et al.l [l996). The 
edge of the GGD 37 flo w is approximate l y loca ted at 
Source S as defined in lHartigan fc Ladal (jl985f h and 
the W2 and S l ocations coincide with a line of soft 
X-ray emission (Pr avdo fc Tsuboll [2005), as well as 
velocity-resol ved [O III] (5007 A) emission in the case 
of Sou rce S ([Hartigan et al.l 1 1 9861 ). iPravdo fc Tsuboil 
(2005) identify W2 and S as separate exciting sources for 
soft X-rays, although they observe significant emission 
from the line connecting the two as well, and the [O III] 
emission has linewidths of ~ 420 km s _1 . It is therefore 
not very surprising that such highly ionized species as 
3+ and Ne 4+ are present; these ions trace the hottest 
gas and provide constraints on the shock energetics. 

Our data support the identification of GGD 37 as a 
protostellar outflow originating at W2. The continuum 
spectral index a is defined as d log(X F\)j d log(X). Cal- 
culated at the position of the Wl-2-3 aperture (Table [T]) 
a ~ 1 when measured from 5.5 to 10 jam and ~ 3 when 
measured from 10 to 35 /im, both results in the range 
typical of Class 0/1 YSOs. Our data reveal spectral lines 
with extremely high excitation potentials, ranging up to 
97 eV in the case of [Ne V]. The ratio of [Ne II] /[Fe II] (26 
/am) which peaks in the vicinity of W2, is 10 - 20 - much 
larger than is seen in other i ntermediate mass YS O re- 
gions (e.g. HH 54 or HH 7- 1 1 lNeufeld~eTal1l2006bf ) . The 
[Ne III]/[Ne II] ratio is ~ 0.2-0.5 in this region, gradu- 
ally falling off to the west before peaking slightly again 
at the edge of the flow, although it should be noted that 
this may be confused by separate excitation mechanisms. 
Ne + ions can be produced in T Tau disks by X-ray ab- 
sorption in the warm disk, and the recombination of Ne 
ions is slow. Ne 2+ has on ly weak charge transf er, and 
Ne+ none at all, with H (jGlassgold et al.l I2007D . This 
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Fig. 1. — Top: near-IR (K'-band) image of GGD 37 (Hodapp 1994) in continuum and H2, with key positions annotated. Bottom: 
Location of the 15" HPBW apertures used to extract individual spectra in GGD 37, overplotted on the [O IV] line map, with corresponding 
positions in the near-IR map. 



TABLE 1 

Surface Brightness by Chemical Species by Aperture 



Species Wavelength (/mi) Module W123 Intensity Jet head Intensity 



O IV] 25.9 LH 5.80 1.60 

NeV] 24.3 LH 0.57 0.34 

Ne III] 15.5 SH 78.9 28.0 

Ar III] 21.8 LH 0.41 



Note. — Surface brightness in units of 10 -6 erg/cm 2 /s/sr of all species at 15" half-power beam width (HPBW) apertures of three regions of 
interest: the W1-W2-W3 emission complex (W123), the extreme NW end of the GGD 37 flow (Jethead), and a selected region devoid of structure 
in the southeastern corner of the map (Continuum Region); see Figure In order to integrate the different modules, we utilized the flux of the [Ne 
II] and the continuum level as a normalization constant; the SH fluxes are thus decreased by a factor of ~ 4.5 from observed values, and the SL2 
spectra are decreased by 1.33 to match continuum levels with SL1. A " — " marking indicates a non-detection or failed linefit. 



is of importance in analyzing [Ne II] in protoplanetary 
disks. In the present case the [Ne III] and [Ne II] emis- 
sion regions are so similar in shape, they are likely to be 
excited the same way, in an extended and unusually pow- 
erful protostellar outflow. The spectra closely resemble 
that of a supernova remnant, although there are several 
differences: the lack of ce rtain lines such as [P II] 32.5 
fim and [Ne III] at 36 fim dNeufeld et al.ll2QQ7D . although 
the latter is probably too faint to be detected in this case. 
Additionally, [O IV] and [Ne V] are not observed in HII 
regions. 

Published models (e.g. iHollenbach fe McKed I1989L 
hereafter HM89) don't seem to account for the observed 
line ratios, however. The first problem comes in measur- 
ing the mechanical luminosity of the flow from cooling 
rates. [O I] 63 fim is modeled as the principal cooling 
line, and [Fe II] and [Si II] are both expected to trace 
[O I]. Predicted intensity ratios of [Fe II] /[O I] and [Si 
II] /[O I] are ~ 0.1 for most values of electron density 
(n e ) and shock speed (v s ), as these lines are also optically 
thin. [O I] intensity is proportional to the total mass flow 
through the shock boundary; using the "W123" aperture 
we calculated an expected intensity of ~ 2 x 10 -4 erg 
s _1 cm -2 sr _1 , and the resulting mass flow rate is 10 9 



no v s for density no. Assuming a fast shock velocity of 
100 km s _1 , and no ~ 2 x 10 4 cm -3 , the mass flow rate 
is 3.3 x 10 -6 M0yr -1 . If instead we perform the same 
calculation with [Si II] (34.8 /im), the mass flow is 5.4 x 
10 _5 M yr _1 , nearly 10 times higher. This is due to the 
ratio of [Si II] /[Fe II] flux > 1 over the entire flow, which 
does not match any conditions in HM89. Although the 
[Si II] and [Fe II] emission are morphologically quite sim- 
ilar, this may reflect differing critical densities, or a high 
Si/Fe abundance compared to solar values. However this 
is more likely because HM89 do not include sufficiently 
fast outflows. 

Additionally HM89 do not predict fluxes of high- 
excitation species such as [Ne V] and [O IV] . The Map- 
pings III Shock Model Library (jAllen et all [2008) pre- 
dicts mid-IR flux ratios, as already noted for the case 
of [Ar III]. The [Ne 11} /[O IV] ratio is - 500 in the W2 
region and 190 at the leading edge of the flow. The ratio 
of [Fe II] (26 /im) to [O IV] is 14.3 and 41.8, respectively. 
[Ne III]/[Ne II] is - 0.21 and 0.15, respectively. [S IV] 
is not detected, but an upper bound to [S IV] /[S III] is 
~ 0.11 and 0.35. The relative lack of [O IV] compared 
to [Ne II] and [Fe II] predicts a large shock velocity of 
500 km/s or more, without a precursor. The lack of [S 
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Fig. 2. Top Left: Map of [O IV] emission from W2 (lower left 
corner) to the edge of the ionized flow (upper right corner). The 
emission peaks at two additional locations in the map: first, at 
the convergence of W1/W2/W3, and second at a point just inside 
the edge of the flow of lower ionization spectral lines. Top Right: 
The same region of sky mapped in [Fe II] 26 /im for comparison. 
Bottom Left: The same plot for [Ne V] emission. Bottom Right: 
[Fe II] 24.5 fim. Although the [Ne V] detection is weak, the spatial 
distribution of the emission suggests that it is real, and spatially 
distinct from the stronger [Fe II] emission. 
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Fig. 3. — Comparison of the peak of bow shock emission in a 
high excitation line, [O IV], and a lower excitation line, [Fe II], in 
terms of normalized flux (vertical axis), as a function of position 
on the flow a xis (rotated 40 degrees), similar to the analysis in 
Raines (2000). The solid curve represents [Fe II]. The dashed line 
represents [O IV]. The two vertical lines line up with the peak 
emission in[0 IV], at the W2 and Jethead apertures. Note that 
the "Jethead" peak is offset from the corresponding [Fe II] peak by 
~ 5", or one full resolution element. 



IV] and relatively high [Ne III] flux seem to indicate a 
slower 200 km s _1 shock, but still without a precursor. 
Mappings III is far too low in density (and therefore in 
recombination rates), and does not fit the observed line 
ratios. 

Modeling the molecular emissio n and the forbidden 
lines simultaneously is difficult. iFlower et al.l ([2003) 
model GGD 37 itself with a 25 km s _1 J-shock with 
a magnetic precursor, which matches the H2 emission 
well but is not nearly fast enough to produce the high- 
excitation lines. In advance of detailed complex modeling 
we can make cruder estimates of shock energetics implied 
by the mid-IR and X-ray emission. 

We can estimate the shock speed from detailed bal- 
ance from excitation by UV/X-ray emission, as a func- 
tion of the ratio of [Ne V]/[0 IV]. The highest tempera- 
ture achieved in the shock is given by the [Ne V] 97 eV 
transition, indicating a temperature of 4 x 10 6 K. We 
can calculate the shock sp eed from the observe d temper- 
ature and luminosity (e.g. IWatson et all 120071 ). Assum- 
ing a solar abundance ratio of O/Ne ~ 8, we derive a 
post-shock temperature of 18,100 K for the observed [O 

IV] /[Ne V] ratio of ~ 10. This indicat es a post-shock 
veloci ty of only 24 km/s. Furthermore, iSchneider et al.l 
(2009) posit that the entire flow is one continuous flow of 
shock-heated material with a velocity in the hundreds of 
km s _1 . The approximate emitting volume is ~ 4 x 10 52 
cm 3 , and the derived n e of emitting plasma is ~ 10 cm -3 
in that region. Thus we find a total shocked mass 0.62 
M for a density of 2 x 10 4 cm -3 . The total momentum 
of the flow is then 124 M km s _1 , assuming a fast shock 
speed of at least 200 km s _1 . This is quite high compared 
to other powerful flows obse rved with Spitzer- IRS (e.g. 
1-2 M km s- 1 for HH 7-ll: lMaret eTaI1l200l h further 
evidence that GGD 37 is unusual. 

If we propose a single emitting region containing all 
three observed ionization states of Ne (Ne + , Ne 2+ , 
Ne 4+ ), supposing the hot postshock gas is uniform in n e 
and temperature T, and the ionization states are in equi- 
librium, we could in theory estimate the physical condi- 
tions. However under these assumptions it is impossible 
for a single region to produce all three lines simultane- 
ously. Thus comparable intensities of all three transitions 
suggests either a photoionizing source or a collisionally 
ionizing source - in either case, likely a powerful and ex- 
tremely hot jet, highlighted by the X-ray emission region 
(Figure H]). 

A detailed discussion of the H2 emission, as well as the 
[Fe II] and other fine structure line maps will be explored 
in subsequent papers. 

4. CONCLUSIONS 

We present the detection of [Ar III], [O IV], and [Ne 

V] fine structure lines, the highest appearance potential 
for optical or infrared lines ever detected in the vicinity 
of a Herbig-Haro flow. Considering this data along with 
0.2-1.0 keV X-ray maps of the flow, we suggest that GGD 
37 is the hottest HH object detected so far. 

The extended nature of the high ionization energy lines 
and soft X-ray emission from the GGD 37 flow suggests 
that it is extraordinarily energetic for a Herbig-Haro ob- 
ject, perhaps an indirect detection of a radiative jet. In 
agreement with models of planar C+J shocks, we de- 
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Fig. 4. — Color contour map of 0.3-1.5 keV X-ray emission of a 1' x 0.8' region in GGD37; overplotted in red unfilled contours is the [O 
IV] Spitzer-IRS emission, overlaid in color intensity contours is the [O IV] emission from W2 (lower left corner) to the edge of the ionized 
flow (upper right corner). 



tect a significantly extended transitional boundary layer. 
This suggests that the radio source W2 is a plausible 
candidate for a high mass embedded protostar driving 
an outflow to the northwest. 
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